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Abstract. Measurement of cosmic microwave background polarization is today a major goal of observational 
cosmology. The level of the signal to measure, however, makes it very sensitive to various systematic effects. In 
the case of Planck, which measures polarization by combining data from various detectors, the beam asymmetry 
can induce a conversion of temperature signals to polarization signals or a polarization mode mixing. In this 
paper, we investigate this effect using realistic simulated beams and propose a first-order method to correct the 
polarization power spectra for the induced systematic effect. 
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1. Introduction 



After the success of COBE (Smoot et al, 1992) and 
WMAP (Bennett et al. [21KJ31 , the Planck mission, to be 
launched by ESA in early 2007, is the third generation 
space mission dedicated to the measurement of the proper- 
' ties of the Cosmic Microwave Background (CMB). About 
20 times more sensitive than WMAP, Planck will observe 
the full sky in the millimeter and sub-millimeter domain in 
nine frequency channels centered around frequencies rang- 
ing from 30 to 70 Ghz (for the Low Frequncy Instrument, 
or LFI) and from 100 to 850 GHz (for the High Frequency 
Instrument or HFI). Of these channels, the seven at low- 
est frequencies — from 30 to 350 GHz, are polarization 
sensitive. 

Temperature anisotropies have been detected by many 
experiments now, the most recent of which detect a se- 
ries of acoustic peaks in the CMB spatial power spec- 
trum (de Bernardis et al, 120001 Hanany et al, [2000, 
Benoit et al, 120031 Hinshaw et al, H0O3), confirm the 
Gaussianity of observable CMB fluctuations (Komatsu 
et al, 120031 though wavelet methods have detected pres- 
ence of non-Gaussianity in WMAP data, Vielva et al, 
2004 ) , and demonstrate the spatial flatness of the Universe 
(Netterfield et al. l2Utj2l Lee et al. l2DuH Benoit et al. l2P051 



Spergel et al. 120031 Spergel et al J2U06|) . This provides com- 
pelling evidence that the primordial perturbations indeed 
have been generated during an inflationary period in the 
very early Universe. The next challenge now is the precise 
measurement of polarization anisotropies and, in particu- 
lar, the detection of the pseudo-scalar part of the polariza- 
tion field (the B modes of CMB polarization) which are 
expected to carry the unambiguous signature of the energy 
scale of inflation and of the potential of the inflationary 
field. The Planck mission will be the first experiment able 
to constrain significantly these B modes over the full sky 
and hence to measure them on very large scales. 



The first detection of CMB polarization at one de- 
gree angular scale of resolution, at a level compatible 
with predictions of the standard cosmological scenario, has 
been announced by Kovac et al (|2002p (see also Leitch et 
al YMffl) . Since then, CBI and CAPMAP have also ob- 
tained significant detection of CMB E-mode polarization 
(Readhead et al. I2TM1 and Barkats et aZ. I2PP5) . More re- 
cently, the iJoomerang (Piacentini et a/ 120051 Montroy et 
a? l20u5l) and WMAP (Kogut et al. [20031 Page et al. l2(J06f 
teams have obtained a measurement of the temperature- 
polarization correlation and E-mode spectrum compatible 
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with cosmological model. No significant constraint on B- 
mode at degree angular scales exist today. 

While the measurement of the temperature and polar- 
ization auto and cross power spectra of the CMB carries 
a wealth of information about cosmological parameters 
and about scenarios for the generation of the seeds for 
structure formation, some near-degeneracies exist which 
require extremely precise measurements. In particular, a 
very precise control of systematic errors is required to con- 
strain parameters which impact these anisotropies and po- 
larization fields at a very low level. 

Many sources of systematic errors are potentially a 
problem for polarization measurements. In particular, the 
shape of the beams of the instrument need to be known 
with extreme precision. In addition, when the measure- 
ments of several detectors are combined to obtain polar- 
ization signals, it is required that the responses of these de- 
tectors be matched precisely in terms of cross-calibration, 
beam shape, spectral response, etc. 

Measurements of Planck telescope beams in the actual 
operation conditions are not to be made on ground. Also, 
there arc no polarized astrophysical sources for the in- 
flight beam calibration. Therefore, one should rely on nu- 
merical simulations of the beams and self-correcting algo- 
rithms of data processing that should allow efficient elimi- 
nation of systematic errors. In polarization measurements, 
a significant systematic error would arise due to elliptical 
shapes of telescope beams which appear, mainly, due to el- 
lipsoidal shape of telescope mirrors introducing astigmatic 
aberrations and other beam imperfections even with the 
otherwise ideal mirror surfaces. 

In this paper, we investigate the impact of beam imper- 
fections on the measurement of polarization power spec- 
tra. We then discuss a method for first-order correction of 
the effect of these imperfections. To illustrate this method, 
we apply it to the case of Planck HFI polarization mea- 
surement. 

The remainder of this paper is organized as follows. In 
section [2j we discuss the issue of beam shape mismatch 
for the detection of CMB polarization. Section [3] is dedi- 
cated to the computation of simulated readouts using the 
realistic beams described in Section [2] and the reconstruc- 
tion of polarized power spectra. In Section |4l we present a 
method to correct for the systematic bias in the B mode 
power spectrum induced by the asymmetry of the beams. 
Finally, Section [5] draws the conclusions. 

2. The beam-mismatch problem in polarization 
measurements 

The HFI polarimeters employ Polarization Sensitive 
Bolometers (PSB) cooled down to the temperature of 
100 mK by a space 3 He - 4 He dilution fridge. These de- 
vices, also used on Boomerang (though at 300 mK), 
are presently the most sensitive operational detectors 
for CMB polarization measurements (Jones et al., 2003, 
Montroy et al, 2003). Each PSB measures the power of 



the CMB field component along one linear direction spec- 
ified by the PSB orientation (Turner et al. [200T|) . 

Ideally, the PSBs are combined in pairs, each pair 
placed at the rear side of respective HFI horn, with the 
two PSB of the pair, a and b, being oriented at 90° of rela- 
tive angle and receiving the radiation from the same point 
on the sky. The ideal polarimeters produce the measured 
signals (readouts): 

s a = -(I + Q cos 2a + [/sin 2a) (1) 

s b = -(I- Q cos 2a - [/sin 2a) (2) 

where /, Q, U are the Stokes parameters of incoming radi- 
ation (see e.g. (Born and Wolf. [TM7|) for the definition of 
Stokes parameters) and a is the angle between the orien- 
tation of the first (a) of two PSBs and the first (a;) of two 
orthogonal axes of the frame chosen for the representation 
of Stokes parameters. The V Stokes parameter does not 
enter Eqs. ([I]), ((2]) since the PSBs are designed to be, ide- 
ally, insensitive to V and, besides, V is extremely small 
for the CMB radiation. 

In practice, the detectors produce beam-integrated 
signals so that equation [T] is modified to (Kraus 1986) 

s a = \ I dfi f dv {IJ + Q a Q + U a U + V a V) (3) 

^ ./beam -/band 

and similarly for s;, where the PSB responses I a {x, v) etc 
are the telescope beam patterns of Stokes parameters com- 
puted in transmitting mode and normalized to unity at 
maximum, functions of both the radiation frequency v 
and the observation point x (the V term is neglected in 
most of the following discussion) . The responses of differ- 
ent polarimeters (a and b) should be adjusted as much as 
possible (both in frequency and in angular pattern on the 
sky) so that, ideally, one should have 



la = lb (4) 

Q a = I a cos 2a, U a = I a sin 2a, (5) 

Qb = — lb cos 2a, Ub = — h sin 2a (6) 

V a = V b = (7) 



where a, similarly to the definition above, is the angle 
of nominal orientation of polarimcter a on the sky with 
respect to the reference axis chosen for the definition of Q 
and U. In this case, with I a — lb, Eq. ^ is reduced to 
the form similar to Eqs. fTJ) and @. 

For simplicity, we approximate the PSB response as 
averaged over the frequency band of the particular chan- 
nel, thus introducing the band-averaged beam patterns 
defined as 

T(x) = I dvl(x,u) (8) 

J band 

and similarly for Q and U (for radiation independent on v 
on the beam width scale, this generates an exact readout). 
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Spin axis 




Fig. 1. (a) The scanning of the polarized detectors pro- 
vides the measurements of intensity of the CMB field com- 
ponents along four directions at each point on the scan 
path, (b) Definition of axis specifying the Stokes parame- 
ters reference frame as seen from the sky. The detector pair 
on the left (e.g. 143-4) measures the Q Stokes parameter, 
while the pair on the right (e.g. 143-2) measures U when 
Q and U are defined with respect to the (x a ^y a ) frame 
(see Appendix lAl for further detail on reference axis). 



Ideally, the beam patterns on the sky I(x) should be 
as close as possible to a perfect Gaussian. Unfortunately, 
design and construction imperfections, telescope aberra- 
tions, and optical misalignment all generate small differ- 
ences in the beam patterns, the impact of which must be 
investigated accurately, especially for very sensitive CMB 
polarization measurements. 

Measuring polarization, i.e. measuring the /, Q and U 
Stokes parameters, indeed involves combining several such 
measurements with different angles a to separate the /, 
Q and U contributions. The Planck HFI detector set-up 
is such that the beams of two horns with complementary 
pairs of PSB oriented at 45° one pair with respect to the 
other follow each other on circular scan paths on the sky 
as shown in Figs. Q] and [2] (see Appendix [X] for further 
detail and notations). Then, in a system where reference 
axes for defining Q and U are along the scan path (x) and 
orthogonal to it (y), the four readouts s a of this set of PSB 
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Fig. 2. Planck focal plane unit (FPU) with polarization 
sensitive bolometers as seen from the sky. Complementary 
pairs of PSB detectors are arranged in two horns following 
each other while scanning the sky so that four detectors 
are in an optimized configuration for polarization mea- 
surement. 

(q = 0°, 45°, 90°, 135°) allow the direct measurement of /, 
Q and U as 



I = Sqo + Sg o = S450 + S1350 

Q = S Q° ~ s 90° 
U = S450 — Si35° . 



(9) 



In practice, when the responses Io° an d /90 are not 
perfectly equal, there is a small residual of / in the esti- 
mate Q of Q: 



Q = Q 



1 



dQ(I Q o - I g0 o)I. 



(10) 



beam 



Similarly, there may be a small leakage of each Stokes 
component into the others. These errors are a source of 
trouble for measuring B mode especially, as they result in 
the leakages of I into E and B (possibly significant) and 
of E into B because I E ^> B on most scales. 

This source of systematic effects for polarization mea- 
surements is not specific to Planck. Any instrument mea- 
suring polarization in a similar way, where signals propor- 
tional to / need to be eliminated from the measurements 
in order to obtain polarization data, may suffer from this. 

The quantitative investigation of the impact of such 
effects requires a realistic estimate of mismatch between 
the companion beams, the simulation of signal data us- 
ing these beams, the reconstruction of maps and of power 
spectra using these simulated data, and the investigation 
of the correction of the effect by data processing methods. 

The computation of the Planck HFI beams is discussed 
in Appendix [Al and in (Yurchenko et al., I2004b|) . Here, we 
will describe their main characteristics relevant for the 
following sections. Because of telescope aberrations, the 
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shape of the intensity beams is essentially Gaussian ellip- 
tic (down to nearly -30 dB) with the major axis around 
10% longer than the minor axis (see Fig. [3]). Thanks to 
the use of PSB, the intensity beams of an orthogonal 
pair of detectors within one horn are very close to each 
other: the difference is at most 0.6%. On the other hand, 
the difference between beams of detectors in two different 
horns can be up to 7%: this difference is mainly due to the 
different orientation of the beam ellipses. As emphasized 
in Appendix [Aj relaxing the assumptions of perfect con- 
ductors and perfect alignment is not expected to strongly 
modify the general shape of the beams. In the next sec- 
tions, we will refer to these computer-simulated beams as 
"realistic beams" . 

3. Effect of beams on polarization power spectra 

The main goal of this section is to study the systematic 
effect induced on the power spectra estimation by realis- 
tic beams described in the previous section, knowing that 
this effect will depend also on the scanning strategy. As 
the Planck mission will scan the sky along large opening 
angle circles, resulting in large parts of the sky where the 
scans are mostly parallel, we have focused the study to 
the observation of a 15° x 15° region of the sky scanned 
only along parallel directions. This restriction does not 
spoil the interest of the study as the small scale distor- 
tion of the beams are expected to affect mainly the small 
angular scales of the power spectra. In addition, other ex- 
periments scanning only a fraction of the sky are affected 
by the similar systematic effects. This restriction also of- 
fers a practical advantage: the computation of the effects 
of tiny beam mismatches on sub-beam scales requires a 
map resolution better than the beam size. For the Planck 
HFI 143 GHz channel, the resolution of about 7 arcmin- 
utes justifies models at sub-arcminute scales. We have thus 
chosen to work on maps of 2048 x 2048 pixels of about 30 
arcseconds each. 

The following paragraphs describe the generation of 
CMB polarization maps from power spectra, and the sim- 
ulation of instrument signals. 

3.1. Generation of CMB polarization maps 

Simulated square maps of CMB intensity and polarization 
are generated using the approximate relation between the 
power spectra in flat (C(k)) and spherical (C;) coordi- 
nates: k 2 C(k) ~ l{l + l)Ci\ l=k (see, for example, White 
et al I1999p . The three maps of T, E and B are then 
computed from three independent realizations of Gaussian 
white noise -D^k), D 2 (k) and D 3 (k) as: 
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Fig. 3. Broad-band power patterns (/ responses) of the 
telescope beams to be superimposed on the sky for polar- 
ization measurements, (a) HFI-143-2a and (b) HFI-143- 
4a, in SC spherical frame on the sky, with the spin axis of 
telescope as a pole (isolevels are shown from the maximum 
down to — 60 dB with a step of — 3dB). 



so that the correlation between the T and E maps is taken 
into account. CVs are the usual spectra describing the 
CMB temperature and polarization. For our simulations, 
we used the cosmological parameters from the WMAP 
best fit model, except that we imposed a tensor to scalar 
ratio of 0.1. The simulated maps include the Gaussian part 
of the gravitational lensing effect of the E mode. 
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3.2. Simulation of instrument readouts 

The readouts must be computed from the /, Q and U 
Stokes parameters. We thus need to convert the E and B 
maps to Q and U using relations (38) in Zaldarriaga & 
Seljak [19971 

a Q (k) = a E (k) cos 20 k - a B (k) sin 20 k (11) 

a u (k) = a E (k) sin 20 k + a B (k) cos 2</> k (12) 

where k x + ik y = fee^ k . The readout from one detector is 
then obtained by convolving its /, Q and U beams with the 
/, Q and U maps from the sky and summing as in Eq. ([3]). 
In the case of the parallel scanning strategy we used, the 
convolution can be easily done, once for all directions of 
observation, by multiplication in Fourier space. Thus, we 
obtain four maps of readout signals, one for each polariza- 
tion channel, so°j S90°j S45° an d Si35°, with polarization 
angles a — 0°, 90°, 45° and 135° with respect to the x- 
axis of the map. With account of established focal plane 
unit (FPU) notation of channels (see Appendix |A|) . they 
correspond, e.g., to the PSB channels 4b, 4a, 2b, and 2a, 
respectively, of two horns HFI-143-4 and HFI-143-2 where 
:r-axis is the ipsc-sxis of spacecraft (SC, see Appendix [A"| 
frame viewed from the sky (Fig. [I})). 

Since the goal of this work is to study only the system- 
atic bias induced on polarization power spectra, we do not 
add any white or low-frequency noise to the signal, nei- 
ther any other systematic effects (Kaplan & Delabrouille 
2002). These other systematic effects will be studied in 
detail in a forthcoming paper. In particular, we assume 
here that the time constant of bolometers, which induce 
an elongation of the beams in the scanning direction, has 
been corrected for. 

3.3. Reconstruction of the power spectra 

The parallel scanning strategy allows us to reconstruct the 
/, Q and U maps from the readout maps using Eqs. ([9]). 
The reconstructed E and B maps can be obtained from 
Q and U using the reciprocal transformation of the equa- 
tions (fTTj) and (fT2"]l . The power spectra are then estimated 
directly from the Fourier transform of the reconstructed 
/, E and B maps, by averaging the a x (k)a Y * (k) in bins 
of width Afc = Al = 20 (with X,Y e {I, E, B}). The re- 
covered power spectra are then corrected for the smooth- 
ing effect due to the beams, which can be approximated 
in Fourier space by a factor exp[— 1(1 + l)er 2 ]. However, 
because of the pixelization of the maps, this approxima- 
tion is not good enough. Instead, we have corrected the 
power spectra using the power spectrum of the intensity 
beam, B(k) = (la^k)! 2 ), where a 7 (k) is the average of 
the intensity beams of the four detectors. This is exact if 
the beams are axially symmetric and identical, and other- 
wise provides a way to symmetrize the beams in Fourier 
space. We have used this correction in all the power spec- 
tra shown hereafter. 



The B mode power spectrum reconstructed by using 
an ideal circular Gaussian beam in both the readout and 
reconstruction computations (assuming Eqs. {5]) and ((6|)) 
is shown on Fig. [4ji. The points shown are the average of 
450 simulations and the error bars represent the disper- 
sion. The relative error is shown in Fig. 2)3, demonstrating 
that the statistical error on the power spectrum recon- 
struction averaged over 450 simulations is less than 2%. 
Finally, Fig. [5] presents the histogram of the bias divided 
by the dispersion, which is well fitted by a Gaussian with 
unit dispersion as expected for the ideal case. Identical 
results are obtained with other power spectra (T, E and 
T-E correlation). 

We can now use this tool to estimate the bias induced 
on the power spectra reconstruction by the realistic beam 
shapes described in section [2] 

3.4. Effects of beams on polarization power spectra 

We apply our algorithm (both the readout simulation and 
C'i reconstruction) using the realistic beam patterns /, Q 
and U presented in section [2j The output power spectra 
shown in Fig. [6] and [Jj are averaged over 450 simulations. 
The temperature power spectrum is perfectly recovered, 
while we can distinguish a small but systematic excess in 
the E power spectrum at I > 2000 and a systematic loss in 
the T-E correlation for I > 1000. The B mode is strongly 
affected after the peak of the lensing signal at I ~ 900, 
with a bias of up to 50% of the signal at I ~ 1500, and 
about 10% around the lensing signal peak at I ~ 1000. 

The spurious B mode may come from leakage of either 
the temperature or the E mode. In order to separate the 
two possible origins, we have done the same simulation us- 
ing the realistic beams /, Q and U when computing both 
the readouts and Cj reconstruction but with no input E 
mode (and no T — E correlation) . The results for T-E, E 
and B power spectra are shown on Fig. [9] (the tempera- 
ture power spectrum is not modified). We observe that the 
spurious B mode is about three times smaller in this case, 
indicating that about 2/3 of the spurious B seen with re- 
alistic input E mode came from E leakage. On the other 
hand, the level of the E mode is much higher than would 
be expected if it came from a mixing from B modes into 
E modes. 

In order to check this, we made a simulation with el- 
liptic Gaussian beams identical for detectors within the 
same horn, but with different ellipse directions for differ- 
ent horns. We have assumed Eqs and [5] to represent 
the ideal polarization sensitivity of the channels, so that 
there is no total intensity leakage into polarization signal 
due to beam difference, and, again, used no input E mode. 
In this situation, the recovered E mode is a small fraction 
of the input B mode signal, i.e. much smaller than the 
recovered E mode in the previous test. This means that 
the E mode recovered in the previous test (using realis- 
tic beams) was not due to a leakage from B modes to E 
modes, but rather from a leakage from T to E. We con- 
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Fig. 4. (a) Input and recovered B mode power spectrum 
with an ideal instrument, i.e. when four identical and axi- 
ally symmetric Gaussian beams are used for both the read- 
out generation and the C; reconstruction. The small peak 
at I ~ 100 is produced by the gravitational waves (the 
tensor to scalar ratio is r = 0.1), while the main pattern 
peaking at I ~ 1000 is due to the lensing effect. The er- 
ror bars are smaller than the thickness of the black solid 
line showing the input model, (b) The relative error be- 
tween the recovered and the initial power spectrum; the 
recovered power spectrum is the average of 450 simulta- 
tions: the statistical error is less than 2%, thus allowing 
the detection, in non ideal cases, of biases higher than 2%. 
Identical figures are obtained for T, E and T — E correla- 
tion power spectra. 

elude that, both the E mode and the spurious B mode 
found with realistic simulated beams when there is no E 
mode in the readout simulation, come from a tempera- 
ture leakage due to the differences in the beam patterns 
between detectors within the same horn, which is up to 
0.9% (Fig.ETJi). 

We thus see that two different effects produce the ob- 
served spurious B mode. First, there is a mixing between 
the two polarization modes, essentially from E to B as 
E ^> B on all scales, due to the beam mismatch between 




Fig. 5. Histogram of the biases divided by the statistical 
dispersion for all multipole bins shown on figure [4j As 
expected, the histogram is well fitted by a Gaussian of unit 
variance, showing that the dispersion on 450 simulations 
gives a good estimate of the errors. 

the two different horns. Second, there is a temperature 
leakage, this time due to the beam mismatch between the 
PSB within the same horn. 

As seen on Fig. [Jj the beam asymmetry affects mainly 
the high-/ part of the power spectra (typically I > 500). 
However, it is not negligible at low I, where the gravita- 
tional waves lie. The Fig. [5] shows the recovered power 
spectra when there is no initial B mode in the simula- 
tion compared to the expected B mode signal from grav- 
itational for various tensor-to-scalar ratios. The leakage 
from T and E mode to B mode power spectrum becomes 
greater than the gravitational wave B mode signal for 
T/S< 10~ 5 . 

3.5. Link with previous work 

Various studies have been done on the systematic effects 
on CMB polarization measurements, in particular the ex- 
haustive work by Hu et al., 120021 (referred to as HHZ 
hereafter). This paper tries to estimate analytically the 
systematic effects on B mode power spectrum, using a 
second order expansion and relating the terms of the ex- 
pansion to beam defects such as, for example, pointing 
error, ellipticity, monopole leakage or calibration. The dif- 
ferent systematic effects are assumed to be described by a 
statistically isotropic field, with a power spectrum of the 
form : 



Ci cx exp(-^ + l)a 2 ), 



(13) 



so that the leakage from T or E to B can be written as 
a convolution between EE or TT and systematic effect 
power spectra (Eqs. 30 and 34 in HHZ). 

In our approach, the /, Q and U maps are given 
through the signals of four detectors (Eqs. [5]). In the quasi- 
ideal case of all Q and U beams defined by Eqs [5] and ® 
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Fig. 6. Input and recovered power spectra of (a) temper- 
ature and (b) T-E correlation signals, using the simulated 
beams of Sec. [5] for the readout simulation and the C; re- 
construction. The recovered power spectra are corrected 
for an average symmetric beam effect by multiplying them 
by the power spectrum of the average beam map. 
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Fig. 7. Input and recovered power spectra of (a) E mode 
and (b) B mode signals, using the simulated beams of 
Sec. [5] for the readout simulation and the C/ reconstruc- 
tion. The recovered power spectra are corrected for an 
average symmetric beam effect by multiplying them by 
the power spectrum of the average beam map. 



but with intensity beams different between the two horns, 
it can be shown that the error on the B mode power spec- 
trum is given by : 



5C? B 



A/i 3 (l) cos 2 (200 sin 2 (200 



B, 



-Cf E {a) (14) 



where A/13 = h — h is the beam difference, Bi is the 
average beam power spectrum and Cf E (a) is the power 
spectrum of the E map convolved with the average beam. 
This form is a particular case of the one given by HHZ, 
and describes completely the systematic effects on the B 
mode power spectrum due to the beam mismatch between 
horns. Note however that there is no mixing between dif- 
ferent I of the power spectra. The major difference with 
HHZ estimate is in the shape of the beam difference power 
spectrum (first factor in Eq. IT4|) which can be fitted by 
Cf 7 oc Z 4 in the interval I e {1, . . . , 3000}, very different 
from HHZ assumption, Eq. [T31 



In the particular case we consider (beam mismatch) 
our approach gives more realistic results, as we use accu- 
rate simulations of the beams. Moreover, the power spec- 
trum of the defect due to beam mismatch we compute 
from these beams is very different to HHZ assumption, 
leading to a different estimate of the size of the effect. 

4. Correction of polarization spectra for 
systematic errors 

We shall now propose a simple way to correct for the spu- 
rious B mode deduced from the observations of the previ- 
ous section. The idea is to assume that temperature and 
E mode maps are recovered well enough to estimate the T 
and E to B mode leakage, if we know the beam patterns. 
We will discuss three cases of C\ correction, depending on 
the knowledge of the beams. In all three cases, the initial 
readouts are generated with the realistic /, Q and U beam 
patterns of Sec. El 
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Fig. 8. Recovered B mode power spectrum in a simu- 
lation with no initial B mode. The theoretical B mode 
power spectra due to primordial gravitational waves are 
also shown for different values of the tensor-to-scalar ratio 
: 0.1, 1CT 3 , 1(T 5 and 1(T 6 from top to bottom. 

4.1. Perfect knowledge of the intensity beam pattern 

In order to have an idea of the ability of the method to 
remove the spurious B mode, we have tested it in the case 
of a perfect knowledge of the intensity beam patterns. 
However, because of the lack of polarized point sources 
with known polarization characteristic, we assumed that 
only the intensity beam patterns / were perfectly mea- 
sured while the Q and U needed for the Ci correction are 
computed using relations ([5|) and ([6]) with the relevant / 
in all the three cases considered. 

The method is as follows. By a quick analysis of the 
data, we would find maps and their corresponding power 
spectra similar to the ones shown in Figs. [6] and [71 Since 
the T and E mode power spectra are recovered with a 
very good approximation, we may assume that the recov- 
ered maps are good as well. Starting with the temperature 
and the E mode maps, assuming no initial B mode and 
using a precise knowledge of the beams, we could then 
simulate the instrument signals. From the previous con- 
sideration, we expect to find, from these simulated signals, 
a spurious B mode polarization coming both from a tem- 
perature leakage and a polarization mode mixing. The B 
mode power spectrum of this simulated signal should be 
an estimate of the spurious B mode. 

The result for the B mode correction, using exact I a 
and lb and assuming relations (|5|) and j6| for the leakage 
estimation, is shown on Fig. 1101 The correction allows us 
to reduce the bias down to less than 1% of the lensing 
signal in the interval 2 < I < 1500. 
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Fig. 9. Spurious generation of E and B modes from tem- 
perature signals using simulated beams of Sec. [2] for the 
readout simulation and the Ci reconstruction, but with no 
initial E mode. 



4.2. Assuming identical beams within the same horn 

In a second case, we supposed that, in order to increase the 
signal to noise ratio, we need to use the signal from both 
detectors within one horn to measure the beam patterns. 
With this method, we would find as beam pattern the 



average of the beams of the two detectors within one horn, 
i.e. the average error on the beams is about 0.5% of the 
beam maximum. 

The result obtained for the B mode correction is pre- 
sented on Fig. QT] This time, there is still some bias left 
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Fig. 10. Recovered B mode power spectrum before (red, 
dashed line) and after (green, dot-and-dashed line) correc- 
tion. The power spectrum is corrected by subtracting the 
estimated leakage (blue dotted line) assuming knowledge 
of the exact beams (see text). Bottom: Difference between 
corrected and initial power spectra. 



Fig. 11. Recovered B mode power spectrum before (red, 
dashed line) and after (green, dot-and-dashed line) correc- 
tion. The power spectrum is corrected by subtracting the 
estimated leakage (blue dotted line) using beams averaged 
within one horn (0.5% error, see text). Bottom: Difference 
between corrected and initial power spectra. 



in the corrected power spectrum, around 3% at I ~ 1000 
and up to 13% for I ~ 2100. 

4.3. Fitting the beams with elliptic Gaussians 

If we have only few point sources or low signal-to-noise 
ratio on signal, we may want to parametrize the beam 
patterns with a function requiring a small number of pa- 
rameters. As an example, we have fitted the four intensity 
beam patterns by elliptic Gaussian. The error of the fit is 
around 2% of the maximum of the beam. 

The result is shown in Fig. [12l together with the dif- 
ference between the corrected and initial B mode power 
spectra. The result is very similar to that of Fig. [11] (us- 
ing horn-averaged beams), though the remaining bias is 
slightly higher. 

This simple method thus seems efficient to recover the 
right height of the lensing effect peak at I ~ 1000. Though 
it is applied here in the case of a simple scanning strategy 
(parallel scans), it should be applicable to any scan strat- 
egy, as soon as the bias estimation is done using the beams 
as precise as possible and the same scanning strategy as 
the real one. 

5. Conclusions 

In this paper, we have shown the effect of asymmetric 
telescope beams on the bolometric measurements of po- 



larization of incoming radiation by considering the case of 
the Planck satellite mission. We have used electromagnetic 
simulation of the optical system (including telescope and 
horns) to compute the main beam shapes of the different 
detectors of Planck. These beams are roughly Gaussian el- 
liptical, with a major axis 10% larger than the minor axis 
and with essentially different orientations of the beam el- 
lipses for the two horns to be combined to measure the 
full set of Stokes parameters, /, Q and U. 

By simulating the scan of a patch of the sky by Planck 
with these realistic, simulated beams, we have estimated 
the bias induced on the E and B mode polarization spec- 
tra due to their asymmetric shapes. We first remark that 
the E mode power spectrum is very well recovered (once 
corrected for an effective symmetric beam), the bias be- 
ing around 0.1% of the signal in the multipole range 
300 < I < 2000, where lies the most interesting part of 
the signal. On the other hand, the B mode is affected 
by a bias around 10% at the peak of the lensing signal 
(I ~ 1000) and increasing for higher I, up to 100% of the 
signal at I ~ 2500. This bias has two origins. First, it is 
produced by the difference of beam patterns of two differ- 
ent horns combined to measure Q and U . This difference 
induces mainly an error on the polarization angle, which 
turns to a mixing of E and B modes. Since, in general, 
E 3> B, we observe finally a leakage from E to B. The 
second origin of the bias is the minor difference of beam 
patterns of two PSB channels with orthogonal polariza- 
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Fig. 12. Recovered B mode power spectrum before (red, 
dashed line) and after (green, dot-and-dashed line) cor- 
rection. The power spectrum is corrected by subtract- 
ing the estimated leakage (blue dotted line) using ellip- 
tic Gaussian beams fitted on the exact beams (2% error, 
see text). Bottom: Difference between corrected and initial 
power spectra. 

tions within the same horn, which induces a temperature 
to polarization leakage. 

Finally, we have proposed a way to correct the B mode 
power spectrum from the above bias in a one-step correc- 
tion which uses the measured T and E maps to compute 
the expected leakage into B when they are observed with 
a model of the instrument's beams. The efficiency of this 
correction depends on the precision of the beam knowl- 
edge: for example, using elliptical Gaussian fits of the ac- 
tual beams allows us to reduce the bias from 10% to 3% 
at the lensing signal peak, / ~ 1000. In all cases, this first 
order correction has been shown to reduce significantly B 
mode contamination. More refined treatments, currently 
being investigated, are expected to be yet more efficient if 
needed. 
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Appendix A: Simulation of the Planck HFI 
telescope beams 

It will not be possible to measure HFI beams on ground. 
The HFI bolometers indeed work only at 100 mK and are 



designed for the thermal load of a few Kelvin environment 
in space. Whereas the instrument can be put in a large 
vacuum tank cooled to 4 K, it is not possible to perform 
far field measurements of the full system, which would 
require placing a source hundreds of meters away from 
the instrument. Hence, responses can only be measured 
at subsystem level (e.g. bolometers + horns) and must be 
associated to a physical model of the telescope to predict 
the beam shape of the complete integrated optical system. 

In this paper, for the investigation of beam mismatch 
effects, we use computer-simulated Planck HFI beams. We 
consider four HFI-143 beams comprising eight PSB chan- 
nels used for the polarization measurements in the band 
centered at the frequency v = 143 GHz. The polarization 
direction of each PSB is specified by the polarization an- 
gle ip on the sky and labeled by the relevant index of the 
channel (0i a , V'lfe e ^ c ) as shown in Fig.QJb). 

In the design of the focal plane unit (FPU), the po- 
larization angles ipi a notifying the channels are measured 
from the upward direction of local meridian of the spher- 
ical frame of spacecraft (SC) having the geometrical spin 
axis of telescope as a pole and counting the angles 
as viewed from telescope, clockwise toward the direc- 
tion of maximum polarization sensitivity of the channel. 
Similarly, we define the polarization angle ip at each ob- 
servation point x in the beam pattern. The direction of 
maximum polarization sensitivity is the major axis of po- 
larization ellipse at point x\ for the angles ipi a notifying 
the channels, x is the beam axis defined as the point of 
maximum intensity / (at this point, the beam field is lin- 
early polarized). 

Following this definition, we consider eight PSB chan- 
nels of the HFI-143 beams which are sensitive to the lin- 
early polarized radiation with polarization angles on the 

Sky lpi a = 2 a = 45°, 1p lb = 1p2b = 135°, -03a = "04a = 0°, 

and 03b = 04b = 90° . The four beams are arranged in two 
pairs (1 and 3, 2 and 4), with two beams of each pair scan- 
ning the sky along the same scan path as shown in Figs.[T] 
and[2j In each pair of beams, the angles 0i Q correspond to 
a full set of four PSB detectors for polarization measure- 
ments with optimized polarimeter configuration (Couchot 
et al. li"999l) . 

The power patterns of two beams tracing the same 
scan path, HFI-143- 2-a/b and HFI-143-4-a/b, are shown 
in Fig. [3] as projected on the sky in the spherical frame SC 
with coordinates tpsc,9sc (vsc = 90° — Osc) where the 
azimuthal angle ipse is counted to the right from the opti- 
cal axis of telescope as viewed from the sky and the polar 
angle Bsc is measured from the upward direction of nomi- 
nal spin axis (the optical axis of the telescope corresponds 
to ipse = 0° and r/ S c = 5°). 

Notice that both the a,b labels of channels and polar- 
ization angles ip are conventionally defined with respect 
to meridians (verticals) of the SC frame viewed from tele- 
scope, as accepted by the FPU design. In the meantime, 
because of the scanning strategy, the reference axes for the 
Stokes parameters on the CMB maps are usually parallels 
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(horizontals) of the SC frame viewed from the sky to the 
telescope. 

To reconcile these definitions, we continue to use the 
polarization angles ip and the established notations a, b 
for the PSB channels. In the same time, for processing 
the readouts according to Eqs. (P)-©, we define beam 
responses /, Q, U, V in SC frame viewed from sky, with 
the first and the second reference axes being the azimuth 
ipse an d the elevation r]sc, respectively (they constitute 
the right-hand frame xy for defining Stokes parameters 
on the sky as viewed from sky to telescope). With these 
definitions, the polarization angle in xy frame is the angle 
a in Eqs. Q, flSJl where a = ip + 90°. 

The beams in Fig. [3] are computed with an extended 
version of the fast physical optics code (Yurchenko et al, 
2001) developed specifically for the efficient simulations 
of the Planck HFI beams. The extended code allows us to 
propagate via the telescope the aperture field of the HFI 
horns mode-by-mode at various frequencies. The aperture 
field is generated by the PSB bolometers considered as 
polarized black-body radiators (in the transmitting mode) 
located at the rear side of the horns. In this way, we obtain 
the band-averaged far-field patterns of Stokes parameter 
responses /, Q, U, V of the broad-band telescope beams 
as produced by the actual corrugated horns (Yurchenko 
et al. I2002a[) rather than by simplified model feeds. 

Rigorous computations of beams require scattering 
matrix simulations of horns (Murphy et al, 2002). In this 
approach, the effective modes of the electric field at the 
horn aperture, E nm , are represented via the canonical TE, 
TM modes £ n j of a cylindrical waveguide as follows 

2M 

E nm (p, V) =22 S nmj£nj(Pi V 1 ) (A.l) 
3=1 

where S nm j is the scattering matrix computed by Murphy 
et al (|2002|) for each horn at various frequencies, n = 
0,1,. ..,N is the azimuthal index and m,j = 1,2, ...,2M 
are the radial indices accounting for both the TE (m, j = 
1,...,M) and TM (m, j = M + 1, ...,2M) modes. 

Recent simulations of the HFI- 143 beams (Yurchenko 
et al, I2004aj) were performed with the scattering matrices 
of size 20 x 20 [M = 10, N = 1) using nine sampling fre- 
quencies spanning the band v = 123 — 163 GHz. Although 
the power patterns of these beams only slightly differ from 
those computed earlier (Yurchenko et al, I2002ap with ma- 
trices 10 x 10 and five sampling frequencies (AP < 0.1 dB 
at P = -3 dB and AP < 1.5 dB at P = -30 dB), 
the effect on the difference between the beams of differ- 
ent polarization and on the fine polarization properties of 
beams is noticeable. This suggests that the latter param- 
eters could be sensitive to other features of the model as 
well. 

In this paper, the HFI-143 beams are computed with 
two essential updates compared to (Yurchenko et al, 
2004a) : the horn design is slightly altered so that the horns 
are now slightly elongated compared to those used earlier, 



and the horn positions are now the final ones, being de- 
fined by the parameter Rc — 1.2 mm that specifies the 
refocus of the horn aperture with respect to the geomet- 
rical focus of telescope for each beam. 

The horn positions were optimized to achieve the best 
resolution (the minimum beam width) of the broadband 
beams (this also maximizes the gain) so that the value 
Rc = 1.2 mm is close to the optimal horn positioning. 
We use updated scattering matrices of size 20 x 20 for 
representing the horn field and nine sampling frequencies 
for spanning the frequency band v = 121 — 165 GHz which 
is characteristic of the updated horns. 

At this stage, we assume smooth telescope mirrors 
with ideal elliptical shape, perfect electrical conductivity 
of their reflective surfaces, and with ideal positioning of 
both the mirrors and horn antennas. The convergence ac- 
curacy of computations was better than 0.1% relative to 
the maximum of the beam intensity pattern I(ip, rf) (for 
comparison, the difference of the broadband power pat- 
tern and the pattern computed at the central frequency is 
about 1%). 

To be perfectly representative of what the actual 
beams may be, one should take into account possible mis- 
alignments of the optics, tilts and deformations of the mir- 
rors, etc. In principle, however, tolerances on the align- 
ment of mirrors and positioning of horns are such that 
the modifications they induce on the beams are supposed 
to be small, and we neglect this last issue for the present 
work. 

To minimize the beam mismatch between the channels 
of orthogonal polarizations, the PSB bolometers of each 
pair of channels, o and 6, are placed in the same cavity 
in the rear side of the respective horn and share the same 
optics — waveguides, filters, horns, telescope. Because of 
this design, the difference of power patterns of orthogonal 
polarizations of the same beam on the sky is really small 
(Fig. IA.11 a), with the peak value for all the broadband 
beams being about dl^aAb = max(/4 a — iAbj/Imax = 0.6% 
(at the central frequency / = 143 GHz, the typical differ- 
ence is 5l4a4b — 0.9%). 

A small difference of this kind arises for two reasons: 
(a) due to minor axial asymmetry of polarized modes that 
appears on the horn aperture when the PSB radiation (in 
the transmitting mode) propagates through the horn (the 
difference varies in sign and magnitude with frequency, 
though being well balanced over the band) and (b) due 
to some difference in the propagation of different polar- 
izations along the same path via the telescope (all the 
differences are computed with the patterns normalized to 
the unit total power of the beams). 

The mismatch of power patterns of different beams is 
about 10 times more significant (Fig. IA.l| b). It depends 
essentially on the location of horns in the focal plane of 
telescope. For the pair of beams HFI-143- 2 and HFI-143-4, 
when superimposed on the sky by spinning the telescope 
until the coincidence of azimuths of beam axes, the peak 
difference of the relative power across the pattern varies 
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Fig. A.l. Relative difference of power patterns (a) I 4a — 
lib of two orthogonal channels of the same beam HFI-143- 
4 (tp — 0° and tp = 90°, respectively) and (b) 7 4a - I 2a 
of the channels HFI-143-4a and HFI-143-2a (-0 = 0° and 
ip = 45°, respectively) when superimposed by spinning 
the telescope about the geometrical spin axis. 



from 7.0% to 8.2% depending on the polarizations being 
compared. Notice that the statistical difference of 5% is 
already rather crucial for the reliable reconstruction of the 
CMB polarization map (Kaplan et al,[2002). 

Fig. shows the patterns of Q and U Stokes param- 
eter responses of the HFI-143-2a and HFI-143-4a beams, 
respectively. The peak values of these parameters are 
Q 2a = 0.6% and U 4a = 1-2% (ideally Q and U should 
be zero in these polarizations). For comparison, the peak 
values of V are 3.7% and 4.2%, Respectively The positive 
and negative values of Q and U (as well as V) are well 
balanced over the beam patterns and the average is very 
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Fig. A. 2. The Q and U Stokes parameters responses: (a) 
Q of the beam HFI-143-2a and (b) U of the beam HFI- 
143-4a (ideally, both Q and U should be zero in the beams 
of these polarizations for the selected reference axes). 



close to zero. It proves that the chosen directions of polar- 
ization of the horn aperture field as found by optimizing 
on-axis beam polarization directions (Yurchenko, 2002b) 
are pretty good, even though the beam patterns are not 
quite symmetrical due to aberrations. 

The analysis of different contributions to non-zero val- 
ues of Q 2 a: U 4a and V shows that V arises mainly be- 
cause of the field propagation via the telescope (VJ nax on 
the horn aperture is only 0.5%). On the contrary, both 
Qia = 0.6% and U 4a — 1-2% given above and the power 
differences between the orthogonal channels of the same 
beams, Sl2a2bj= SliaAb — 0.6%, are essentially due to the 
horn effects (U = 0.4% and 5I a b = 0.8% on the horn aper- 
ture). The telescope contribution, though non-additive, is 



C. Rosset et al.: Beam mismatch in CMB polarization measurements 



13 



still important, as the propagation of the axially sym- 
metric quasi-Gaussian source field shows (in this case, 
Sha2b = 0.9% and £Zi 46 — 1-0% in the beams on the 
sky, being zero in the source field). 

Finally, in the cross-beam power differences Sl2a4f3, 
both the horn and the telescope effects are significant (e.g., 
^^2a4/3 depends, to some extent, on polarizations being 
compared), although the telescope effect dominates (for 
the quasi-Gaussian source field, 8l2a4/3 varies from 5.7% 
to 6.7% in a way consistent with the variations in the 
beams from the actual corrugated horns of respective po- 
larizations) . 
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